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CONNECTION COEFFICIENTS FOR CLASSICAL 
ORTHOGONAL POLYNOMIALS OF SEVERAL VARIABLES 

PLAMEN ILIEV AND YUAN XU 


Abstract. Connection coefficients between different orthonormal bases sat¬ 
isfy two discrete orthogonal relations themselves. For classical orthogonal 
polynomials whose weights are invariant under the action of the symmetric 
group, connection coefficients between a basis consisting of products of hyper¬ 
geometric functions and another basis obtained from the first one by applying 
a permutation are studied. For the Jacobi polynomials on the simplex, it is 
shown that the connection coefficients can be expressed in terms of Tratnik’s 
multivariable Racah polynomials and their weights. This gives, in particular, a 
new interpretation of the hidden duality between the variables and the degree 
indices of the Racah polynomials, which lies at the heart of their bispectral 
properties. These techniques also lead to explicit formulas for connection co¬ 
efficients of Hahn and Krawtchouk polynomials of several variables, as well as 
for orthogonal polynomials on balls and spheres. 


1. Introduction 

For a positive measure p defined on that satisfies some mild assumptions, the 
space of orthogonal polynomials of degree n in d variables with respect to the 
inner product 

{f,9)=[ f{x)g{x)dp[x) 

has dimension dimV^ = where 0 P G if P is a polynomial of degree 

n and (P, Q) = 0 for all polynomials Q of degree less than n. The elements of 
a basis {P^, : v G Nq, \v\ = n} for may not be orthogonal among themselves. 
A basis whose elements are orthogonal to each other is called orthogonal, that is, 
(Pjy, P^) = 0 for V ^ pL, and it is called orthonormal if, in addition, (Pj^, P^) = 1. 

If := {Py : V G Nq, |u| = n} and Qn := {Qi/ : u G Nq, \v\ = n} are two 
orthogonal bases of then we can express one in terms of the other; for example, 

(1.1) Qi , — ^ |u| — n, G R. 

|/i|=ri 

We call the coefficients connection coefficients of in terms of P„. Through¬ 
out the paper we adopt the convention of using P^ = P,y/||Pi/|| to denote the 
orthonormal polynomial when P^ is an orthogonal polynomial, and we denote by 
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the normalized connection coefficients between two orthonormal bases. Thus, 
(HU can be rewritten as 


Qv — ^ ^ 

|/j|=n 


If both {Pf^} and {Q^} are orthonormal bases, the matrix C := is orthogonal 

(0 p. 67]) and, therefore, satisfies 

( 1 - 2 ) ^ ^ and ^ ^ 

\uj\—n \u)\—n 

We are interested in identifying the connection coefficients for classical orthogo¬ 
nal polynomials of both continuous and discrete variables. Especially interesting 
is the case when can be expressed explicitly in terms of discrete orthogonal 
polynomials and their weights. 

Our starting point is the Jacobi polynomials of two variables that are orthogonal 
with respect to the weight function — x — y)'^ on the triangle = {(x, y) : 

X > 0,y > 0,1 — x — 2/ > 0}. For this weight function, one orthonormal basis 
of consists of polynomials Pj{x,y) := P^!^P{x,y) for 0 < j < n given in 
terms of the classical Jacobi polynomials. Another orthonormal basis consists of 
Qj{x,y) := P^"Jj°‘j{y,l — x — y), obtained from Pj by applying the permutation 
(123) on (a,/?, 7 ) and (x^y,l — x — y) simultaneously. It was shown in [1] that the 
connection coefficients in this case are expressed in terms of the Racah polynomials 
(see Section[3l) by taking a limit of the corresponding results for Hahn polynomials, 
whereas the proof of the latter result in [3] uses specific irreducible representations 
of the symmetric group. 

We start with an elementary proof of the result in [3] for two variables and then 
we consider the connection coefficients for the Jacobi polynomials of arbitrary num¬ 
ber of variables on the simplex for bases generated by elements of the symmetric 
group. The picture in several variables quickly becomes much more complicated. 
The constructions require new ingredients and the results lead to interesting phe¬ 
nomena. For instance, in the case of three variables, the connection coefficients can 
be expressed by the Racah polynomials, the Racah polynomials of two variables 
defined in or by the sum of products of such polynomials, depending on what 
permutation is applied in obtaining the second basis. Moreover, in the case of the 
Racah polynomials of two variables, the dual orthogonality provides a natural in¬ 
terpretation of a subtle duality between the degree indices and the variables, which 
lies at the heart of the bispectral property studied in [5]. Indeed, in our framework, 
the duality amounts to transposing the orthogonal matrix C and comparing the 
orthogonality relations in (11.21) . All this is trivial for the one variable Racah poly¬ 
nomials, for which the symmetry between the variable and the degree index follows 
immediately from the explicit hypergeometric representation. 

From the results about the Jacobi polynomials on the simplex, connection coeffi¬ 
cients for several other families of classical orthogonal polynomials can be derived. 
First of all, using the Jacobi polynomials on the simplex as generating functions of 
the Hahn polynomials [HIIH], the connection coefficients for the Hahn polynomials 
of several variables and their permutations can be derived from those for the Jacobi 
polynomials on the simplex. Secondly, taking an appropriate limit ([nmH]) , we 
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can derive explicit formulas for the connection coefficients of the Krawtchouk poly¬ 
nomials of several variables, which can be expressed by Krawtchouk polynomials 
with different parameters. Finally, it is known that the orthogonal polynomials on 
the unit ball and the unit sphere and those on the simplex are related, which can 
be used to derive connection coefficients for orthogonal bases on these domains. In 
particular, this leads to explicit formulas for the connection coefficients for orthog¬ 
onal bases given in Cartesian coordinates and in polar coordinates, respectively, on 
the unit ball. 

It should be pointed out that this work is connected to several other areas. The 
explicit formulas of the connection coefficients for the Hahn polynomials of two 
variables in [3] and those for the Jacobi polynomials on the simplex in [4] were 
re-proved in [14] from Lie and quantum algebra representations; see also [15] and 
the references therein along this line. Classical multivariable orthogonal polyno¬ 
mials have also appeared recently in a large body of work around the quantum 
singular oscillator model and, more generally, quantum systems, which appear to 
share many common points with our work, see [HIT] and reference therein. In 
the language of representation theory and mathematical physics, our connection 
coefficients share common trait with the Clebsch-Gordan coefficients and the 3nj 
symbols. While we adopt a more traditional special function approach, we believe 
that it would be interesting to find a Lie-theoretic interpretation of the new results 
concerning the multivariable Racah polynomials in the present paper. Note that for 
the large family of multivariable Krawtchouk polynomials introduced in [9], such 
Lie interpretation together with their spectral properties were obtained in [10] . For 
a recent introduction to these polynomials and numerous probabilistic applications, 
see [5] and the references therein. 

The paper is organized as follows. In the next section we recall the Jacobi poly¬ 
nomials on the simplex, define the connection coefficients for these polynomials and 
their permutations, and study their general properties. Since the explicit formulas 
for these coefficients are given in terms of the Racah polynomials of one or several 
variables, we collect the results for the Racah polynomials in Section [31 including 
several new results clarifying the dual Racah orthogonality. Explicit formulas of 
the connection coefficients in dimension 2 are given in Section [J] which serves as 
a basis for further study, and Section [5] contains a detailed study of all connection 
coefficients in dimension 3. In Section jh] we present the general results in arbitrary 
dimension and, in particular, we give the explicit formula for the cyclic permutation 
in terms of the Racah polynomials of several variables. The connection coefficients 
for the Hahn and Krawtchouk polynomials are discussed in Sections |7| and jS] re¬ 
spectively. Finally, the connection coefficients for orthogonal polynomials on the 
unit ball and the unit sphere are studied in Section |9j 


2. Connection coefficients for Jacobi polynomials on the simplex 

The Jacobi weight function on the simplex := {a; € > 0, |a;| < 1}, 

where |a;| := xi -I- ... -I- Xd, is defined by 


a; € 


( 2 . 1 ) 


W,,(x) := a:”^ • • • x'^‘‘ (1 - |a;|)'"'^+i 
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where k G with ki > —l,...,Kd+i > —1. The Jacobi polynomials of d 

variables are orthogonal polynomials with respect to the inner product 

(2-2) (/.g)vE. := -I-d+i / fix)9{x)W4x)dx. 

The normalizing factor in front of the integral is chosen so that (1,1)^ = 1. Let 
denote the space of orthogonal polynomials of degree n with respect to 
Wk,- a family of orthogonal polynomials with respect to W^, can be given explicitly 
in terms of the Jacobi polynomials. For a, fd > —1, the Jacobi polynomial 
satisfies 


(2.3) 

(2.4) 




—n,n + a + /3 + l 1 — x 


a + 1 


(a + 1)„ f 1 + X 


n\ 


2 Pi 


—n, —n — / 3 _ X — 1 
a + 1 ’ X + 1 


To describe our orthogonal polynomials on the simplex, we will need the following 
notations that will be used throughout this paper: For y = (yi,... ,yd) G and 
1 < J < d, we define 


(2.5) Yj :={yi,...,yj) and := {yj,... ,yd), 

and also define yo := 0 and := 0. It follows that y^^ = y^ = y, and 

\yj\ = yi + ---+yj, |y^l = % + •••+ yd, and |yo| = |y^+^| = o. 

For K = (ki, ..., Kd+i), we have := [kj, ..., Kd+i) for 1 < j < c? + 1 and adopt 
this notation also for other Greek letters. For u G Ng and k G we define 

(2.6) aj ■.= aj{K,v) ■=\k^^^\ + 2\u^^^\+ d —j, ^ < j < d. 

Notice that ad = Kd+i since |zz‘^+^| = 0 by definition. The standard basis on the 
simplex is given in the following proposition [5l p. 150]. 


Proposition 2.1. For k G with Ki > —1, u G Ng and x G R'^, define 


(2.7) 


p;(x) :=n(l-|xi-i|)‘^^Pi; 

i=i 


2xi 


1 - |Xj_l 


- 1 


where Oj = aj{K, u) is defined in (EH). The polynomials in {Pfi : juj = n} form an 
orthogonal basis afV^fW^.) with A^{k) = (P”,Py) given by 


( 2 . 8 ) 


Ai,{k) := 


(|k| + d + l) 2 |i/| 


n 

i=i 


{kj + Qj + l)2vi (kj + l )^-(Oj + l)y. 


(Kj + flj + 1 ) 




We shall call the polynomials the Jacobi polynomials on the simplex. 

Let 5(i+i be the symmetric group consisting of all permutations of d+1 symbols. 
For T G Sd+i and y G R'^+^ we define 

•“ (yT(i)i... I yT(d+i)). 

To capture the symmetry of the simplex, we make the following definition. 

Definition 2.2. For t G Sd+i and x G we define rx as the first d components 
o/r(xi,... ,Xd,Xd+i), where Xd+i := 1 - jxj; i.e. rx := (Xi-(i),... ,Xr{d)) € T‘^. 
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Evidently the weight function satisfies Wtk{tx) = Wk(x). Furthermore, it 
is easy to see that the set of polynomials {P^'^{tx) : |p| = n} is also an orthogonal 
basis of V^(VE„). Moreover, the square of the norm of is A^{tk). Indeed, 

using the fact that WK(r“^ 2 /) = Wmiu), 

[ [K'^iTx)f Wk{x)dx = [ [P^'^iy)f Wrn{y)dy = A^{tk). 

Jrd JTd 

Since {P^ : |z/| = n} and {PJ’”(t{-}) : \v\ = n} are both bases of V^(Wk), we can 
write one basis in terms of the other. 

Definition 2.3. Let r G Sd+i and p,/ r G Ng with \v\ = |/r| = n. The connection 
coefficients are defined by 

( 2 . 9 ) P^irx) = Y. cIJk)P;:{x). 

\n\=n 

Let Pffix) = Pffix). Then {Pfi : \v\ = n} is an orthonormal basis of 

V^. We can write (12.91) as 

(2.10) P^'^iTx) = Y ’^ith := J 

\li\=n V ^ 

where cJ^(k) are the normalized connection coefficients. 

The connection coefficients satisfy the following useful relation. 

Proposition 2.4. Let ti and T 2 be two elements of Sd+i- Then, for yL,v G Ng and 

W\ = \d\ = n, 

(2.11) Y 

\uj\—n 

Moreover, formula (12.111) holds also for the normalized connection coefficients. 

Proof By dM]), we have Pf'^'^{T 2 x) = Applying n on both 

sides we get 

Pf^^^^{TiT2x) = Y 

\oj\—n 

= Y cD%(n«;) Y 

\uj\—n \^\—n 

= Y c);%(riK)c))(^(K)p;^(x), 

\fiL\—n \uj\—n 

from which the identity (12.111) follows from the dehnition of cffjfiffi). The proof of 
(12.111) for the normalized coefficients follows along the same lines by working with 
formula (12.101) . □ 

As an immediate corollary we obtain the following useful relations for the nor¬ 
malized connection coefficients. 

Proposition 2.5. For p,G Ng with |z/| = |/i| = n we have 

( 2 . 12 ) Y Y 

|u;|=n \u)\—n 
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and 

(2.13) c;_/(k) = 

Directly from the orthogonality, it follows that 

We are interested in finding an explicit expression for It turns out that 

these coefficients can be regarded as polynomials, which satisfy, by ( 12 . 121 ) . a dis¬ 
crete orthogonality relation. In fact, these discrete orthogonal polynomials can be 
expressed by the Racah polynomials. 


3. Racah polynomials 


Since the explicit formula of the connection coefficient cj ^{k) involves the Racah 
polynomials, we collect what we need for these polynomials in this section. 

The classical Racah polynomials in one-variable are defined by 

i?„(A(x)) := Rn{X{x);a,P,'y,S) = 4 F 3 

for n = 0,1,..., and A(a;) = x{x -I- 7 -I- <5 -I- 1) and one of the bottom parameters 
is —N, see [H p. 190]. They satisfy the orthogonality relation 

N 

x—0 

where r„ is the square of the norm of Rn and wr is the weight function given by 


—n, n -|- a -|- ,5 -|- 1 , —a;, x + ^ + 5 + 1 . 
a-l-l,/3-l-(5-|-l,7-|-l 


(3.1) Wfl;(a;) :=w(a;;a,/3,7,(5) 

_ (7 + + l)a;((7 + + 3)/2)a;(Q; -|- l)x{!d + + l)a:(7 + 
a;!((7 -1-5-1- l)/2)a;(7 -|- 5 — a -I- l)a;(7 — /? + l)x(5 -I- l)x 

For d > 2, the Racah polynomials of d variables are defined on the simplex 
Vat = {a; € Ng : 0 < xi < a ;2 < ■ • • < Xd < N} and they are orthogonal with 
respect to the weight function 


(3.2) 


wr{x-I3,N) := P 
1=0 


(/3j + l Pj^Xj + l—Xj ^Pj + l)xj + l+Xj 

{Xj +1 — Xj)\{f5j -|- l)xj+i+Xj 


n 

i=i 


m + 2 )/ 2 ) 

{fXj/2')xj 


where Xg = 0 and x^+i = N. A basis for V„(ixr) is given explicitly by {R,y : u £ 
Ng, \v\ = n}, where 


Ri.(x;/ 3 , iV) y[( 2 |i>'j_i| + fij — MujiWj-i] + Pj+i + Xjj^i)^^ (l^^i-il ~ ^3+1)1- 

1=1 

p + P3+1 - ^0 - 1 , l*^i-il - Xj, + Pj + a^i. O 

I3j — P3,\vj_i\+Xjj^i,\vj_i\—Xjj^i 

These polynomials are mutually orthogonal with respect to ixr, that is, 

^ R^(x; /3, iV)R/,(x; /3, A^)n;R(x; /3, N) = 5^,,^ [ri,(/3, N)f , 
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where N) is the norm of Ri/(-; /3, N). Using formula (2.4) in [17], one can show 
that 


(3.3) [r.(/3,iV)]2 


{Pd+l)N+\i^\{ — N)\u\{ — ^ — /3o)|i/|(2|p| + Pd +1 — f3o)N-\i^\ 
Wo + l)Ar 


d 

X ]^Pfc!(/ 3 fc+i — /lfc)j/j,( 2 |i/fc_i| + / 3 fe — /?o)j/fc(|i^fc| + Wk-i\ + Pk+i — Po — ^)vk- 
fc=i 


Remark 3.1. When d = 1, the formulas above define the same polynomials, up to 
some unessential factors. More precisely, if we set 

(3.4) a = -N-l, /3 = ^ 2 -/ 3 o -1 + A^, 7 =/3i -/3 o - 1, ^ =/3o, 

it is not hard to see that weights in (13.11) and (13.21) differ by a factor independent 
of X and the 4 U 3 functions defining the polynomials coincide. In particular, if we 
consider the orthonormal polynomials we have R„(x) = (—l)”.R„(A(a;)). 


Next we formulate an important duality relation which lies at the heart of the 
bispectral property of the Racah polynomials studied in | 8 ]. Let us define dual 
indices p, variables a), and parameters P by 


(3.5) 


Xj = N-\ud+i-j\ for j = l,...,d, 

Vj = Xd+ 2 -j - Xd+i-j for j = 1 ,..., d, 

Po = Po, 


Pj = Po- Pd+2-j - 2 iV + 1 for j = 1 ,..., d + 1 . 


Proposition 3.2. The map {x, P,v^N) —> (x, P,i', N) is an involution. Moreover, 
the Racah polynomials satisfy the following duality relation 

/g Qx _ R,,{x;P,N) _ 

(_iV)|,|(-iV - /3o)|h UUM +1 - WU 

_ _ Ri^{x;P,N) _ 

" (-iv)|,|(-iv - Po)ioi uUCW+i - Wh, ■ 


Proof. The duality can be obtained from the one established in [ 8 ] , by applying the 
involutions Xj —>■ —Pj — Xj to the duality relations (4.1) and Theorem 4.4 there. 
Alternatively, we give a direct proof as follows. Applying the Whipple identity 


(3.7) {UU{V)m{WU 4 F 3 


(-m,X,Y,Z \ 

V U,V,W ’V 


{1-V + Z- m)m 


y.{l-W + Z 


7R)m(bt)m 4 T 3 


f -m, U - X,U -Y,Z 
\l-V + Z -m,l-W + Z -m,U 



withm = Vj, X = Vj+2\uj-i\+Pj+i-PQ-l, Y = \uj-i\-Xj, Z = \vj-i\+Pj + Xj, 
U = Wj-il - Xj+i, V = 2\uj_i \ + Pj - Pq,W = \ + Pj+i + Xj+i, we see that 

the jth term in the product defining the Racah polynomials can be rewritten as 


(l^i-i| Xj+i)ujW Wj\FPo + Xj)ujW Pj Pj+i+Xj Xj+i),^^ 

^ p (-^p-Xj+i + Xj,-Xj+i - \vj\ - Pj+i + /3o + 1, Wj-i\ + Pj + Xj ^ 

^ Xj+i,l-\Vj\+Po+Xj,l- Vj + Pj - Pj+I +Xj - Xj+I )' 

Plugging the dual variables (13.51) in the last formula, one can show that the 4^3 
term coincides with the 4 T 3 term above with j replaced by d+1 — j. Thus all 4 T 3 
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terms in Ri/(a;;/3,iV)/Rt,(x;/3,TV) cancel and after we simplify and rearrange the 
remaining products, we obtain equation ()3.6|) . □ 


(3.8) 


Using formulas ()3.2I1 . (I3.3|l and (I3.5|l one can check that 

ri,(/3, N)'^ wr{x] P, N) {l3d+i)2N0d+i)2N 


(-TV)|,|(-iV - /3 o)|h nti(/3fc+i - 


[iV! iPo + l)iv] 


which is a constant independent of u. 

Combining Proposition [3^ with the last equation we obtain the following corol¬ 
lary. 

Corollary 3.3. The orthonormal Racah polynomials satisfy the duality relation 

(3.9) \/ wr{x] /3, N) p,N) = \Jwr{x\ P, N) Ri>(i; P, N). 

Remark 3.4. Following Tratnik [H], let us define new variables, indices and param¬ 
eters by 

x'j = N - Xd+i-j for j = 1,..., d, 

(3.10) i'j=i2d+i-j for j = l,...,d, 

/?' = -2N - Pd+i-j for j = 0, 1,... , d + 1. 

Then it is easy to check that 

(3.11) R,{x-p,N) = K,{x'-p',N), 
where R' is defined by 

d 

Kix-,P,N) := l[{2y+^\+Pd+i-l3jU 
i=i 

X {\iy^+'^\ - N - - Xj-i)^^{\v>^+^\ - N + xj-i),,^ 

^-Vj,Vj + 2\u3+'^ \ + Pd+i - Pj-i - 1, \ - N + Xj,\ - N - Pj- Xj . 

2fy^+i| + Pd+i - Pj ,- TV - /?,_i - x,_i, -N + a:,_i 


X 4^3 


Since the Racah weights wr{x' , /?', N) and wr{x, /?, N) differ by a factor independent 
of a:, we see that the last formula defines a second family of Racah polynomials of 
d variables. In particular, if we combine this with Corollarv l3.31 we obtain 


(3.12) 

fy wr{x; /3, N) R^(x; j3,N) = \^ 

Iwr{x';P',N) %,f. 

where 


II 

for j = 1 ,..., d. 

(3.13) 

lyj — X Xj 

II 



~ + l —Po — 
^(,+1 =-2TV -/3o. 

1 for j = 0,..., d. 
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4. Jacobi polynomials of two variables 


In this section we consider the Jacobi polynomials on the triangle, or simplex 
when d = 2 . In this case, the basis of V^(Wk) consists of Pj, for Pi + 1/2 = n. 
For later use, however, we will write v = {n — j,j), so that the polynomials in (12.71) 
become 

(4.1) P„%.^.(xi,X2) = 

for 0 < j < n. We use the usual notation for permutations to denote the elements 
of the symmetric group S3, 


53 = {(1),(12),(13),(23),(123),(132)}. 

Let r be an element of S3 and x = {xi,X2) G By Definition 12.21 tx = 

ixr{i),XT(2))^ where X3 = 1 — xi — X2- To simplify the notation, we define 

'^’ 3 —(^ 1 1 ^2 ) ■ ^n—j,ji^T{l)y^T{2))- 

Accordingly, we denote the connection coefficients by cj^(/c,n) and (12.9|) becomes 

n 

(4.2) P^i';-^(xi,X2) = cl^{K,n)P^_^ ,^{xi,X2). 

m—0 


Under permutations of the basis Pn-j,j ; there are essentially three distinct fam¬ 
ilies of orthogonal polynomials. The first one consists of P^-j j given above, the 
second one consists of 


(4.3) Pii"]’nxi,X2) = Pi%+"^+"^+^’"^^(2:r2 - 1)(1 - X2)^Pf 


2xi 

1 — a ;2 



for 0 < j < n, and the third one consists of 
(4.4) 

Pi-]’-ixi,X2) = -2X1 -2a;2)(xi +I2)^pf 


2^2 

Xl + X2 



for 0 < j < n. Indeed, there are essentially these three families since it is easy to 
see that 

^’i!_y’”(a;i,a;2) = {-ly P^_jj{xi,X2), 

(4.5) Pn-jj''ixi,X2) = i-iyPi^Yj ixi,X2), 

P^-i^{xi,X2) = (-l)^Pn’;(^G^ 2 ). 

Since (123) = (13) (12) and (132) = (12) (13), the last two identities can also be 
derived from and P^Yj'j^ which explains our notation tk, that is, writing 

the group action on the left side. 


Lemma 4 . 1 . For 0 < j < n, 

(4.6) Fulfil, X2) ^(_i)n-.-(^2 + l)n-,-(^3 + l), ^ (-n) (n+H + 2)^ 

ml{K2 + l)m 

^ f-m,m +K2 +K3 + l,-j,j + Ki +K3 + I \ , 

-n,K3 + l,n+\K\+2 
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Proof. By (I2.3|l with P^’^\—t) = (—and (12.4L it is easy to see that 

={-iy 


3(i2);'^n ^ -( (^2 + l)n-j{nz + l)j 

in- j)\j\ 

—n + j,n+ j + |k| + 2 
K2 + 1 ’ 


X 2 F 1 


X 2 ) 2F1 


- 3 , -3 - n.1 

K3 + 1 


\X2 


Writing the product of the two 2 ^"i as a single sum and rearranging the terms by 
using 


_ i-3)ni{-ni)i j-j - K.i)m-i _ j-j - Ki)mi-K3 - m)i 

m!(l-m + j)i {ks + l)m-i (^3 + l)m(l - w + j + ki)* ’ 

we see that the product of the two 2 ^"i is equal to 


^ i- 3 )m{-j - Ki)m /'-n + j,n+j+\K\+2,-m,-K3-m \ 

2—' mliK3 + l)m ^ \ K2 + l,l-m + j,l-m + Ki+j ’ y 2- 

We can think of each term in the above sum as a meromorphic function of j, which 
is analytic at the nonnegative integers. This 4 T 3 can be rewritten by using the 
iterated Whipple identity [5J (4.5)] 


iu)miv)miwU4F3(^ 

= {1 - X - m)m,(l -v + y- m)m(l -v + z- m)ra 

^ p ( —m,w — x,u — x,l — v — m 
^ ^ \1 — X — m, 1 — V + y — m, 1 — V + z — J 

with X = —j, y = m + K 2 + K 3 + 1, z = j + Ki + K 3 + 1, u = —n, v = K 3 + 1, 
and w = n + \k\ + 2, which leads to the 4 F 3 in (14.61) with a constant that can be 
combined with the constant in front of the original 4 F 3 . Putting these together 
proves (ITO . □ 


We are now ready to derive explicit formulas for the connection coefficients. 
Proposition 4.2. For 0 < j < n, the connection coefficients for r = (12) satisfy 


(4.7) cf^^{^,n)=DlM4F3 

where the constant is given by 

(4.8) DIM =(-l) 


—m, m + K 2 + K 3 + 1 , —j, j + Ki + K 3 + 1 


—n, K 3 + 1 , n + |k| + 2 


;1 


_/ 1 ^rt+m ( n)j {k2 + l)ri-j (^3 + 1) j 


jl(K2 + l)m 

(n + |k| + 2 )^ 


i^2 T ^3 4“ 2.m + 2) n—m iK2 + K3 + m+ 1 ) 

m 


Proof. Setting xi = 1 shows that Pn-m m(l: ^ 2 ) is equal to a multiple of x™, where 
the constant is determined by and the leading coefficient of 

pJk 3 .k 2 )^ More precisely, 

u n\ DK ^ _ («2 + K 3 + 2m + 2 )„_m(K 2 + K 3 + TO + l)m 


{n — m)\m\ 


2 • 
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Hence, the right hand side of (j4.2ll when xi = 1 is a polynomial of degree n in X 2 , 
so is the right hand side, by ()4.6|) . of (I4.2|l . Comparing the coefficients of in the 
two sides proves the stated result. □ 


The 4^3 in the lemma is a Racah polynomial of degree j in m or a Racah 
polynomial of degree m in j. More precisely, let us state the connection coefficients 
in terms of the orthonormal Racah polynomials. 


Corollary 4 . 3 . For 0 < j < n, the connection coejficients for r = ( 12 ) satisfy 
(^) = ai)Rj{X{m); ai), X{m) = m{m + K2 + K3 + 1 ) 

where ui = (—n — 1 , n + ki + K3 + 1, K3, K2), and 

= (-1)"+™+-? y/wR{j;a2)RmiXij);cr2), A(j) = j{j + + K3 + 1 ) 

where (J2 = {—n — 1,n + «;2 + K3 + 1,K3 ,ki). 

Proposition 4 . 4 . For 0 < j < n, let Then 


(4.10) 








Proof. By (14.5L only the case (13) needs a proof. Directly from the definition and 
using the fact that ( -^^2 — l) = (_i)jpj'^ 2 ,Ki)^_^ _ follows that 

piyy;{i,-x2) = (-i)^(23)„pi5;(i,x2), 

where {23)K,g{K) = g{{23)K) and the subscript k indicates that the transposition is 
acting on k. Moreover, gH) shows that 


= (-1)’"P«_^,„(1,X2) 

is invariant under the action of (23). Consequently, setting xi = 1 and replacing 
X 2 by —X 2 in (lO) and applying (23) on k, we obtain 


(-l)^PSd = E cS((23)«)(-l)-P„"_^,^(1,012), 

m—0 


which shows that c^^ ((23 )k') = (—1)'^+™c^-^,^^(k) by the definition of c^^Ik). □ 


Since (12)(13) = (132), the relation (12.111) can be used to derive a summation 
formula between the Racah polynomials. Recall that the weight function for Racah 
polynomials is given in (ED. We define a related weight function w* by 


w(a;;a,/3,7,7) 


w*{x;a,l3,^,S) 


(7 + l)x 
(<5 + l)a; 
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Proposition 4.5. Let u{x; K,n) := w*(x;—n — l,n + \k\ — ki + 1, ki). Then 

n 

(4.11) (-l)’"M(m; At, 11 ) 4^3 

m—0 

XiFs 

_ / ^sn+k+i (^2 + l)k(K2 + 1)£ (ki + K 3 + 2 )„ 

(k 2 + l)ra + l)fe(K3 + 1)^ 

— k, k -t- Kl K2 1, —^ 1^2 4“ 4” 1_ 

—n, K2 4-1, n 4 -|k| 4-2 ’ J' 

Proof. To derive this identity, we use Proposition 12.41 with ti = (12) and T 2 = (13), 
so that (j 2 . 11 l) becomes 



—m, m 4- Kl 4- K 3 4- 1, —k, fc 4- ki 4- K 2 4-1 
—n, Kl 4- l,n 4- |k| 4- 2 

—m, m 4- Kl 4- K 3 4- 1, —£ 4- K 2 4- K 3 4- 1 
—n, K 3 4 - 1 , n 4 - |k| 4 - 2 


= V 4^i2((12)K,n)c|^^^^(K,n). 


m—0 


By ()4.7I1 and (I4.10|) . all three connection coefficients in this identity are explicitly 
given in terms of 4^3 series, which are the three 4 T 3 functions in the statement. 
We still need to check the coefficients in front of 4 F 3 . For the sum in the right hand 
side, these coefficients are (—1)^+'"Z1^^((12)(23 )k)Z1))j ^(k) which can be split into 
a product of two terms, the first term is independent of m and mostly cancels out 
with the terms in Z1^^((23)k) in the left hand side, whereas the second term is 
equal to 

_(n 4- |k| 4- 2)m _n!(K 2 4- l)n-m 

(ki 4- K 3 4- 2m 4- 2)„_m(Ki 4- K 3 4- m 4- l)m (n — m)!m! 

which can be written as u(m; K,n)(K 2 + l)n/(Ki 4- K 3 4- 2)„ using the relations 
= (-l)™(a)„/(l - a- n)m and ( 0 ) 2 ™ = We omit the 

details. □ 


One particular interesting case is ki = K3, for which the weight function becomes 

u{x; K, n) = w{x; —n — 1, n 4- |k| — ki 4- 1, K3, ki) 

and the two 4 F 3 series in the summation of (14.1111 are Racah polynomials of degree k 
and i, respectively, in the variable m associated with this weigh function. Without 
the (—1)'" in the summation, the left hand side of (14.111) would be zero when k ^ i. 


5. Jacobi polynomials of three variables 

In this section we consider the case c? = 3, which is complex enough to give us 
an idea of the complications in high dimensions and also simple enough that we 
can have a clear picture for all permutations. 

In this case, the polynomials in (12.71) are of the form 

(5.1) P;(a;i,ai2,a^3) = /^ir’"^^(2a:i-l)(l-a;i)"=p(f’"=) 

2X3 




I — Xl — X2 
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for |p| = Pi + P 2 + P 3 = R-, where 

0,1 = ai(p, k) = K2 + ^3 + ^4 + 2 p2 + 2^3 + 2, 

02 = 02 (p, k) = K3 + K4 + 2 p 3 + 1 . 

Let r be an element of 54 and x = {xi,X 2 ,X 3 ) G By Definition 12.21 tx = 
(a;T-(i),a:T-( 2 ),a;T-( 3 )), where 2:4 = I — X 1 —X 2 —X 3 . To simplify the notation, we define 

’ (xi,X 2 ,X 3 ) . (^r ( 1 ) j ( 2 ) 5 ( 3 ) ) ■ 

There are a total of 24 elements in 54 . Under permutations of we end up with 
12 different bases, which are given by 

(5.2) (1), (12), (13), (14), (23), (24), (123), (124), (132), (142), (13)(24), (14)(23), 


and bases derived from other permutations can be written in terms of them. Indeed, 


(5.3) 


Pj234);«:(^) ^ 

pa234);«(^^ = (-1)‘'3P^123);«(2;), 


pa2)(34);K(^) ^ 

P,(243);«(^) ^ (_l)‘^3pj24);«(^^^ 

pa243);«:(^) = (-1)’'’^ {x), 


pa324);«(^) ^ (_l)>^3pa3)(24);,^(^)^ P,(1342);«; ^ (-l)'^3pjl32);«^^^^ 

pa423);«(^) = (-1)"^P,(14)(23);.^^)^ P^1432);.(^) ^ (-1)"3P^142);«^^)^ 

The polynomials P” possess inherited structures of two variables, which can be 
utilized to derive formulas for connection coefficients from those for orthogonal 
polynomials of two variables. We start with Cv^fi {k). 


Proposition 5.1. For \v\ = l/rj = n, 

(5.4) c[,^^)(k) - P 3 ) with k = (ki, K 2 , K 3 + K 4 + 2 p 3 + 1) 

= <^>^ 3 .M 3 -D”“;((^(ki,K 2 ,K 3 ) 

p f ~h' 2 , At2 + «^2 + K3 + 1, —P2, P2 + + K3 + 1 , \ 

^ ^ —n + P3, K3 + 1, n + P3 + |k| + 3 ’ ) 


where cl. 2^^2 
(5.5) 


(k, n - 
^(34) 


- P 3 ) are the connection coefficients in (|4.7p . Furthermore, 

(«) = = (-i)'^^c(; 2 )(«). 


Proof. With K3 = 02(k, p) = K3 + K4 + 2p3 + 1 and k = (ki, K2,K3), we can write 

= P 1 . 2 (^ 1 ^^ 2)(1 - xi - X2rp^:^’^^^ (r^ - 0 ’ 

yl - a;i - 2:2 / 

where P^^ = P))_j^ 3 _j ,3 ,^2 is the orthogonal polynomial in two variables given by 

631). Furthermore, the transposition (12) gives 

Pii2);-(x) = Pii,if (a;i,X2)(l - XI - X2)'^='Pi3«-'=^) (r^ - 0 ’ 

yi - xi - X 2 y 

which has a common factor, the one indexed by P3, with P” and, as a consequence, 
it is easy to see that cl^p, {k) = 0 ii V 3 ^ If P3 = ^3, then the connection 
coefficients clearly reduce to those of c\}^}i^{ji,n — P3) for orthogonal polynomials 
in two variables, which proves (Ol) . 
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The other two cases follow readily from the first two identities in ()5.3I1 . 
second one also follows from Proposition 12.41 


The 

□ 


The polynomial P” has another structure of two variables, for which we start 
with a simple observation. 


Lemma 5.2. For r in the subgroup {(1), (23), (24), (34), (234), (243)}, the eonnee- 
tion coefficients satisfy 

(5.6) i^i), k := {k 2 , K 4 ) 

where ^1) connection coefficients for orthogonal polynomials on 

the triangle for the variables [y 2 , 2/3,1 — 2/2 ~ 2/3)- 

Proof Setting 2/2 = and 2/3 = so that 1 - 2/2 - 2/3 = it is 

easy to see that we can rewrite as 

(5.7) Pffix) = Pi“-«^)(2a.i - 1)(1 - a:i)"-‘^^P;,,,3(2/2,2/3), 

where is the orthogonal polynomial in two variables given by 

631). Together with the fact that ai(u, k) = |k| + 2^2 + 2r'3 + 2 is invariant under 
permutations in the subgroup, it follows that r in the subgroup only acts on P^ 2 ,y 3 - 
As a result, we conclude that = 0 if 7^ 2^1. If ui = 211, then (15.71) shows 

that the connection coefficients can be derived from those that express Pf^%^{y 2 , Vs) 
in terms of P”^ ^^3(2/2, Vs) in two variables. □ 

Remark 5.3. In the right hand side of (15.61) . r is acting on (k 2, K3, K3) as a subset 
of (ki, K 2, K3, K4), not on the positions of the elements in (k 2, K3, K4). For example, 
(23)(k2, K3, K 4 ) = (k3, K2 , K4)- 

As a consequence, we can derive from (14.71) and (14.101) the following corollary: 


Corollary 5.4. For \v\ = \pL\ = n, 

(5.8) K3, K4) 


X 4P3 


/~/43i M3 + ^3 + ^4 + Ij ~U3, U 3 + K2 + K 4 + 

y —n + ui, ft:4 + 1, n — Pi + |k| — Ki + 2 



and 


(5.9) 




Furthermore, the identities in (15.31) show that 

= and («) = (-l)^^c(2^)(«). 


In the cases that we deal with so far, the connection coefficients are Racah 
polynomials just like in the case of two variables. This will change in our next case, 
for which it is more convenient to give the formula for the normalized connection 
coefficients. 


Proposition 5 . 5 . For \v\ = n and \fj,\ = n, 

(5.10) c}^^^)(k) = (-l)”+'"3 ^ R(^ 3 _^ 3 ) (x; /3, n), 

where (3 = (ki, /ti + K4 + 1, ki + K 3 + K4 + 2, |k| +3) and x = (P3, P2 + vs)- Moreover, 

(5.11) («:) =c<^i^,3)((132)«:) = (-1)"+^V«^r( 2/;n) R(.3,..) (?/;/?*, n), 
where j 3 * = (k 3, K3 + ^4 + 1, K 2 + K 3 + K 4 F 2 , |k| + 3) and y = (m 3,M2 + Ms)- 
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Proof. Since (123) = (12)(23), we use Proposition [23J which shows that 


\uj\—n 


Since, by (I5.8L {{\2)k) contains and, by (15.4L c2'^J{k) contains <5^3./ij, 

it follows that 

((12 )k)c£^.^_^^(k) with w* = (pi, n - Pi - ^3, ^3). 

Hence, it follows from (15.81) and (15.41) that 


f.( 123 )('^'| _ jjn-ui 


(ki, K3, K2, K3 + K4 + 2^3 + 1) 


X 4^3 


X 4^3 


/'-H3, H3 + K3 + K4 + 1, -P3, P3 + Ki + K4 + 

—n + Pi, K4 + 1, n — Pi + |k| — K 2 + 2 



/-M2, M2 + |k| - Ki+ 2^3 + 2, -n + Pi + /13, n - Pi + M3 + |«:| - «^2 + 
y —n + M3; ^3 4” ^4 4” 2^3 4“ 2,74 4“ M3 4“ |7c| 4“ 3 



Comparing with the explicit formula, it is easy to see that the two 4P3 functions 
are exactly those two that appear in R^(a;;/3,n) with x and /3 as specified in the 
statement. The explicit formula of R^(x;/?,n) also contains a (—X2 )ii 3 (|k| — K 2 4- 
2^2 4- 2)^3. As a result it is sufficient to show that 


-P"3,M3(^0 ^3, ^4)A>n-((i%3,it3(Ki, K2, K3 + K4 + 2fl3 + 1) ^ ^ 

{-n + Pi)ii 3 (|K| - K2 4- n - Pi + 2)^3\/Ai.(K(123))u;R(a;,/3,n) ^ 

where x = (P3, n— pi). By (12.121) and the orthogonality of the Racah polynomials, it 
is sufficient to show that A is independent of x. This involves tedious computation. 
For the record let us write down the intermediate steps. 


A,((123)k) = 


(/t2 4 - 1 )[/i(k 3 4 - 1)1/2 (k 4 4 " 1)1/3 (ki 4 " 1)1/3 
pi!p2!p3!(|k| — K2 4 - 2n — 2 pi + 2 )(ki 4 - m 4 - 2p3 + 1) 
(|k| — K 2 4 - 3)271-1/1 (^1 4 - K4 4 - 2)„ — 1/1+1/3 _ 

(|ai| 4 “ 3)211—1/1(1^1 ^2 4 “ 2)7i_f/j^^j/3 (avi 744 1)1/3 


where P2 = n — pi — P 3 and = means that the equality holds up to a multiple 
constant that is independent of p. Using this formula and the explicit formula of 
wr{x; /3, n), we can verify that 


a/ Ai/((123)7t)wR(x, /3, n) 


(k 2 4- l)i/i {k4 4- 1)1/3 (^3 4- 1)11—1/1—1/3 

vi\v3\{n - Pi - P3)! 


The formula of £>”3,;(3K«;i; «3, «;4)£’((-((i-113,112 «^2, ^3 + ^4 + 2 m 3 + 1) can be 

derived from (14.81) . which can be used to show, together with the above formula, 
that A is independent of p. We omit the details. 

Finally, since (132) = (123)“^, the formula for follows from the one for 

by equation (I2.13|) . □ 


Recall that (k)],/,^ is an orthogonal matrix. Formula (15.101) naturally con¬ 

nects the second orthogonality in (12.121) to the orthogonality of the Racah polyno¬ 
mials. The hrst orthogonality relation in (12.121) is more subtle and can be related 
again to the Racah polynomials via the duality established in Proposition [321 We 
formulate this precisely in arbitrary dimension in Theorem 16.31 
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Explicit formulas for connection coefficients corresponding to other permutations 
can be derived from those of and 

Proposition 5.6. For \v\ = n and \fj,\ = n, 

(5.12) 

((34)«). 

Proof. Since (34)(123)(34) = (124) and (34)(132)(34) = (142), the two identities 
in the first line follow from H2.11I) and c^^)I^(k) = (— 1 )^ 3 ^^^^_ The other identities 
follow from the relations in (15.31) . □ 

All these connection coefficients are given in terms of the Racah polynomials of 
two variables. The remaining permutations are (13), (14), (13)(24), (14)(23), (134), 
(143), (1324), (1423). By (15.31) . we only need to deal with the first four. Let us 
consider, for example, (13). 

Proposition 5.7. For \v\ = n and \fj,\ = n, 

n-V3 

(5.13) = (-1)'"" ^^.A‘(^)-^!-2(AW;o')R(M3.A‘2)((i"3,^ + i^3);/3,n), 

i=o 

where = {-lYyJwR{t, a)yJwR{{v^,i + 1 ^ 3 ), P,n) with 

cr = {—n + U 3 — 1 , n + U 3 + |k| — K 3 + 2 , ki + K 4 + 2^3 + 1 , K 3 ), 

/3 = (ki, Ki + K 4 + 1 , Ki + K 3 + K 4 + 2 , |k| + 3). 

Proof. Since (13) = (123)(12), it follows from p.lip that 
cli2)((l23)K)c^JJ)(/t) 

|ci;|=n 

uii+u:2=n — i'3 

which can be written as a sum over i = uj 2 , since n — uji = uj 2 + vs- D 

Remark 5.8. Using (12.1111 and explicit formulas for connection coefficients, we can 
also derive various relations for Racah polynomials of one and two variables in the 
spirit of Proposition 14.51 For example, using (12) (123) = (23), we end up with a 
formula that expresses a Racah polynomial of one variable as a sum of a product 
of Racah polynomials of two variables and Racah polynomials of one variable. 

Remark 5.9. Note that (15.1311 resembles the formula for the Wigner 9j symbols 
as a sum over triple products of 6 j symbols. It would be interesting to explore 
this connection in more depth and to relate the results in the present paper to Lie 
theory along the lines of the works [mill]. 

6 . Jacobi polynomials of d variables 

Most of the explicit formulas in the previous section can be extended to arbitrary 
dimension. For j € N, let denote the subgroup of Sd+i fixing the points of 

{1, 2,..., j}. Then similarly to Lemma [5.21 we see that the connection coefficients 
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for T G can be computed from the connection coefficients of Jacobi polynomials 
in d — j variables. 

Proposition 6.1. For T G we have 


n — 




— \vj\) are the connection coefficients for Jacobi poly- 


( 6 . 1 ) 

where c. 

nomials on the simplex T‘^~^ in the variables yj+i,yj+ 2 , ■ ■ ■ lUd with parameters 
= (kj+ 1 , Kj+ 2, • ■ •, Kd+i). Moreover, formula (16.1|) holds also for the normal- 


K- 


i+1 - 


ized connection coefficients. 


Proof. The proof of (|6.ip is similar to the proof of Lemma 15.21 The fact that ()6.1I1 
holds also for the normalized connection coefficients follows by using the explicit 
formula for the norms of the Jacobi polynomials (12.8p and the connection between 
given in equation (12.101) . □ 


For fc < d + 1 we think of Sk as the subgroup of Sd+i which acts on the first k 
elements while keeping {fe + l,fc + 2,...,d+l} fixed. Similarly to ProDOsition l5.ll 
we can reduce the computation of for r G S'fc C iSd+i to the computation of 

the connection coefficients in lower dimension as follows. 


Proposition 6.2. For k < d and t & Sk we have 

(6.2) |iz'=+^|), 

where cf,,_ ^^{k,n — are the connection coefficients for Jacobi polynomials 

on the simplex in the variables xi,X 2 ,.. ■ ,Xk with parameters 

k = (ki,K2 , + 2|iz''+^| +d - k). 

Moreover, formula (16.21) holds also for the normalized connection coefficients. 

For the cyclic permutation we obtain the following explicit formula in terms of 
the Racah polynomials extending Corollary 14.31 and Proposition [5^ 

Theorem 6.3. For iz,/x G Ng with \v\ = n and |/x| = n, 

(6.3) ' = (-l)’"+‘^-i y/ wr(£>; ffi n) R/i(£>; /3, n) 

where Ry{y \ ffin) are the orthonormal Racah polynomials of d — 1 variables with 
indices jl = (/i^, p-d-i, ■ ■ ■, 1 x 2 ), variables v = (|iz‘^|, ..., and parameters 

(3j = Ki + | + j for j = 0,1,..., d. We also have 

(6.4) ^n) ffi n), 

where fi,n) are the orthonormal Racah polynomials of d — 1 variables with 

indices v = {vi, ... ,Vd-i), variables fi = (|/i.2MP2li ■ • ■ ; |Pd-il) and parameters 
/3o = Ki, ffi = -2n-d + j for j = l,...,d. 

Proof. The proof of (16.311 can be done by induction similarly to the proof of Propo¬ 
sition [5^ We outline the main steps below. Suppose that the statement is true in 
any dimension less or equal to d and we want to establish it in dimension d J- 1. 
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We can decompose the cyclic permutation as (12.. .d + 1) = tiT 2 G Sd+ 2 , where 
Ti = (12 ... d) and T2 = (d, d + 1). By Proposition 12.41 we have 

(6.5) = E 

\Lo\—n 


Applying Proposition 16.II with T2 = (d, d + 1) G 5^+2 


Using the last formula and the induction hypothesis (or equivalently, Corollary 14.31 
and equations (13.41) 1 we see that 

(6.6) cl^+')((12 ... d)n) = (-1)^^ 

X \/wR(uc;+i;/3',Ud + Ud+i) (ud+i;/?', Ud + Ud+i), 


where /3' = (ki, ki + Kd +2 + 1, ki + Kd+i + Kd +2 + 2). Applying now Proposition [521 
with Ti = (12 ... d) G S'd C Sd +2 we deduce that 

(6-7) - /id+i), 

where k = {ki, K 2 , ■ ■ ■, Kd, Kd+i + Kd +2 + ‘2^d+i + !)■ From the last two equations 
it is clear that the sum in (16.51) can have only one nonzero term corresponding to 
w with coordinates 


UJj = 


Ud + Vd+l — dd+1 
.dd+1 


if j < rf - 1, 
if j = d, 
if j = d + 1. 


In the rest of the proof we fix a; as above and therefore formula (j6.5ll becomes 

(6.8) ... d)K)^^l- '^\K). 

By the induction hypothesis we can rewrite formula (16.71) as 

(6-9) ^wr{v"] P”, n - ^Md+i) P”, n - ^id+i), 

where = {^d, fJ-d-i, ■ ■ ■ , ^2), - kd+i,- /id+i, ■. ■, \v^\ - ^J.d+l), 

/3q = Ki, /3" = Ki + \K‘^+‘^-^\+2fid+i + i+j,j = 1,2,... ,d. Plugging ([63 and (|63 

in (16.8|) . we see that '(i}+"'^~^^\K) contains a product of d 4F3-series and it is not 
hard to show, using the explicit formulas above, that they coincide with the product 
of the 4F3-series defining the d dimensional Racah polynomial on the right-hand 
side in (16.31) for the cyclic permutation (12 .. .d-|-l). Moreover, if WR{i>, I3,n)\d^d+i 
denotes the wr{i>, (3, n) with d replaced by d -I- 1, one can check that 

_ WR{i>,/3,n)\d^d+i _ 

WR{iyd+i]/3', Vd + Vd+i)u)R{iy"]/3",n - ^id+i) 

can be rewritten as 

(Urf -|- Vd+l ~ Md+l)- (^1 4“ l)i^d+i'd+i — dd+l _ 

{Kd+l + Kd+2 + ^)vd+i'd+l+IJ.d+l (^1 + ^d+1 + ^d+2 + ^)i'd+iyd+l+IJ.d+l 

up to a positive factor independent of v. Finally, one can show that the square of 
the norm (u^+i; /?', Ud -|- Ud+i)|p coincides with the inverse of the right-hand 
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side in the last formula, up to a factor independent of p, and therefore 
WR(i>,/3,n)|dH^d+i ||R^^^^(pd+i;/3',Pd + Pd+i)|P 
WR{i'd+i;P',i'd + J^d+i)wR{i^";^'',n- fid+i) 
is independent of p. This combined with the second orthogonality relation in (j2.12l) 
proves formula (16.311 for (12... d+ 1), completing the induction. Finally, equa¬ 
tion (16.41) follows from (j6.3ll and Corollary |3l3l □ 

Remark 6.4. Combining Theorem 16.31 with Bemark 13.41 we also have 

(6.10) /3', n) R(,, (^'; jS', n) 

where R(,, (/i'; /?', n) is the second family of orthonormal Racah polynomials of d — 1 
variables with indices v' = {ud-i,Vd- 2 , • ■ ■, pi), variables = (|/x‘^|, ..., |/.t^|) 

and parameters /3' = \+ j for j = 0,1,..., d — 1, /3^ = —2n — ki. 

Remark 6.5. Note that = {—lY'^P^{x). Using Propositions 16.1116.21 

and Theorem [Ql we see that for every cyclic permutation of the form r = (j,j + 
1,... ,k — l,k), the corresponding connection coefficients can be written explicitly 
in terms of appropriate Racah polynomials in fe — j variables when k < d and d — j 
variables when fc = d-|- 1. For instance, for the adjacent transpositions {j,j + 1) we 
obtain the following explicit formulas: 

^wrYj+i■ (A(pj+i); a) 

where a = {—Vj — Vj+i — l, j, |K:^+^|-|-2|i/'^+^|-|-d—j —1, Kj) 

for j = l,...,d—1 , and 

For more complicated permutations, we can apply Proposition 12.41 and express the 
coefficients as sums of products of Racah polynomials. Furthermore, Remark 15.81 
applies in d dimensional setting. 


7. Hahn polynomials of several variables 
Let be a positive integer. We use homogeneous coordinates of , where 
Z^+i :={a:ae<+Ma| = N}. 

For K G with «:i>—l,l<z<d-|-l, the Hahn weight function is defined by 

(k -|- 1)x 


(7.1) 


H«;,Ar(a;) = 


X G Zl+1 


xl 


'N ’ 


where 1 = The Hahn polynomials are orthogonal with respect 

to the inner product 

(7.2) E /(«)9(«)h..~(o). 




The Jacobi polynomials on the simplex can be used as a generating function for 
one orthogonal basis of the Hahn polynomials. Let Pl^ be defined as in (12.711 and 
let 


pt - n T’, 




1=1 


(i)=n 

1=1 


{ajY,v) + l)^ 

vY 
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Definition 7.1. Let k G with Ki > —1 and N G N. For v G Ng, |u| < N, 

define the Hahn polynomials Hj^ (a; N) for a G by 

(7.3) Pu.N{y) = \y \^—= X! — 

where y = {y',yd+i) G K‘^+^ 


These polynomials are given explicitly in terms of the classical Hahn polynomials 

n, n + a + 6 + 1, —x_ 


Qn{x;a,b,N) := 3 F 2 


1 


a + 1,-7V 

Proposition 7.2. For x G and v G Ng, \v\ < N, 


n = 0,1,..., iV. 


(7.4) H^(a;;K, A^) ^ + \xj-i\ + \iy^+^\), 

X Q^.{xj-,Kj,aj,N - |xj_i| - y+^\), 


where Oj = aj{K, v) are defined in (ED). The polynomials in {Hi,(x; k,N) : \v\ = n} 
form an orthogonal basis o/V)((HK,Ar) and := (H,^(-; k, iV), Hi,(-; k, iV))^ ^ is 
given by, setting '■= |«;| + d + 1, 


(—l)l'^l(AK)jV+|i/| TT (Kj + ttj + l)2l/j(Kj + l)t'3^j! 

(~A^)|[/| (AK)Ar(AK)2|;y| {Hj + + P}uj{a,j + l)iyj 


The explicit formula for H[^(-; k, N) and its norm were stated in |12j by inductive 
formulas. The generating function was later identified as the Jacobi polynomials 
on the simplex and the basis was given explicitly; see [IS [II] for further references. 

A useful observation is that B,^{k, N){p'^)'^ and A;,(/t) = differ by a 

constant that depends only on |u| (see [T8]b 


Corollary 7.3. For v G Ng, 


(7.5) 


B,,{k,N) = 


(_l)kl ( 1^1 d + l)jv+|;y| Ai^(k) 
(-Af)l^l (|/t| + d+ l)Ar (p;5)2 ' 


Remark 7.4. The definition of P” in this paper differs from the one used in 
the constant pj); more precisely, Plf{x)/p'^ is defined as P” in [18]. 


m by 


Let r be a permutation in Sd+i- It is easy to see that {H;y(r{-}, tk, N) : |u| = n} 
is also an orthogonal basis of V)((HK,Ar)- Hence, we can consider the connection 
coefficients, h)) j^(k), of the Hahn polynomials 

(7.6) H,,(rQ;;rK,fV) = ^ k, N). 

|/i|=n 


By the orthogonality, it follows immediately that 

1 


K,M = 


B^(k, A7) 


(Hi-(t{-}; tk, N), H^(-; k, N))„ 
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Let k, N). Then : |/i| = n} is an 

orthonormal basis of V^(HK,Ar). We can write (17.61) as 


Hy(Ta;TK, A^) = ^ 
\fj,\=n 


A^) with h^ ^(«;):= 


B^(k, A^) 


U,fl \ 




[k). 


These connection coefficients are closely related to those of the Jacobi polynomials 
on the simplex, as shown in the following result. 

Theorem 7.5. For z/, /x € Nq and r € 5'd+i, 

(7-7) KA'^) = = 

Pv 

Proof. It is known [151 Cor. 3.5] that, for a € and x S we have 

_ (k + l)a H^(a; K, N) Plfjx) 

{\K\+d+l)N B^(K,Af) 


where X = (x, 1 — |x|). Since (tX)“ = X'^ using the above identity and (17.31) 
with y = tX , we obtain 


prirx) 

Pi'" 


= —rH^(a;TK,X)(TX)“ 


= Y —rH,,(a;T/c,X) 


H” l)Ta 

(|/s:| d l)iv 


E 


H^(t ^a\K,N) Pff[x) 
B^{k,N) pf 


Since the summation is over all a such that |a| = N, permuting the order of the 
summation by r and using (ra)! = a!, it follows that 


Pi'" 


(^ + l)a u / Ar\ H^(a; K, X) T’^(x) 

|K| + d+l)jv a! B^(z< 


^ AT V ’ / 


M<N ’ 

P^i^) 


X) pf 


\t^\<N 

Puix) 

T. 

l„l<« 


Consequently, (1771) follows directly from the definition of cI^^{k). Finally, using 
dLSl), it is easy to verify the identity for hJ_^(K). □ 


Corollary 7.6. The connection coefficients h^^(K) of the Hahn polynomials are 
independent of N. 

In particular, for the cyclic permutation, these coefficients can be given in terms 
of Racah polynomials, as shown in Theorem 16.31 
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8. Krawtchouk polynomials of several variables 

For p G with 0<pi<l, l<i<d and |p| < 1, the Krawtchouk weight 
function of d variables is defined by 


(8.1) K,.„(x):=iV!n ^'-_Mj_, W < JV, 

and the Krawtchouk polynomials of several variables are discrete orthogonal poly¬ 
nomials with respect to the inner product 


( 8 - 2 ) if^9)K,,N ■= fix)9ix)Kp^Nix). 

xen^:\x\<N 

Using the notation (12.51) for p, u and x we can define a family Ki.(-; p, N) of the 
Krawtchouk polynomials as follows (see m)- 

Proposition 8.1. Let p with 0 < p^ < 1 and |p| < 1. For v € Nq, \v\ < N , 
and X G define 


(8.3) 


K^{x;p,N) := 


(-")n 


nviv 


|Xi-l| 


D + 1 I 


X 2F1 


'3^ 


1 - |P,-ll 




P 3 


The polynomials {Ky(-;p, iV) : \v\ = n} form an orthogonal basis of V^{Kp^N) and 
C,y{p, N) := (Kj,(-; p, TV), p, TV))^^ ^ is given by 


C.(p,iV) 


(-l)IH 

(-^)h P? 


where Vd+i = 0. 


Note that when d= 1, we obtain the classical Krawtchouk polynomials K„ {x\p, N) 
of one variable, defined by 

K„(a:;p, A^) := 2 ^"! ^ n = 0,1,..., A^. 

These polynomials satisfy the following discrete orthogonality relation on {0,1,..., A^} 


(8.4) 


-/V x(A \N — 

N\ X; K„(x;p, iV)K„(x;p, 


n\{N — n)!(l — p)" 
N\p^ 


The Krawtchouk polynomials in (j8.3|) can be obtained as limits of the Hahn 
polynomials in (17.41) . More precisely, setting k = t(p, 1 — |p|), we have (see [TT]) 


(8.5) lim H^(x; t{p, 1 - |p|), A^) = (-l)l‘^l N). 

t^oo (1 _ 

We note, however, that the polynomials in (18.31) differ by an unessential factor from 
the ones in m- With the current normalization, the Krawtchouk polynomials 
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polynomials satisfy a simple duality relation. Indeed, if we define dual indices v, 
variables x, and parameters p by 


( 8 . 6 ) 


Xj — ^d-\-l—j 

~ ^ Pd+i-j-(i - IpI) 


for j = l,...,d, 
for j = l,...,d, 

for j = l,...,d, 


then the following proposition holds (see [SI Section 5.4]). 


Proposition 8.2. The map (x, p, p) —>■ {x,v,p) is an involution. Moreover, the 
Krawtchouk polynomials satisfy the following duality relation 

(8.7) K{x; p, N) = Kp{x; p, N). 

From the explicit formula in (18.61) for the dual parameters and by using partial 
fractions it is easy to see that jpj = jpj. Combining this with the formulas for the 
norm and the weight of the Krawtchouk polynomials, one check that 

(8.8) C,(p,iV)Kp,^(x) = (l-|p|)^. 

As an easy consequence of the last formula and the duality in Proposition 18.21 we 
obtain the following corollary. 


Corollary 8.3. The orthonormal Krawtchouk polynomials satisfy 
(8.9) ^Kp^Nix) Kj,(a:; p, N) = ^Kp,iv(i) Ks(i; p, N). 


Let T be a permutation in Sd+i- We denote by rp the first d components 
of t(p, 1 — jpj). Then {Ki,(t{-}; rp,iV) : \v\ = n} is also an orthogonal basis 
of V)((Kp,^). Hence, we can consider the connection coefficients, k))j^(K), of the 
Krawtchouk polynomials 

(8.10) K,,(rQ;;rp, iV) = ><i:,p(p)Kp(a;p,iV). 


By the orthogonality, it follows immediately that 

Let kp{-;p,N) = [Cp{p, p, N). Then {Kp(-;p, A^) : jpj = n} is an 

orthonormal basis of V)((Kp_jv) and we have 


(8.11) Ki,(rQ!;Tp, A^) = ^ 

|At|=rj 


kI).^(p)K^(a:;P,fV), k^ p(p) := 


Arp,N) 


ip)- 


These connection coefficients can be derived from those for the Hahn polynomials. 
Indeed, the limit relation (18.51) leads to the following lemma. 


Lemma 8.4. For r € Sd+i and \i/\ = jpj = n we have 
(8.12) K {p) = lim K {t{p, 1 - IpD). 

In particular, the connection coefficients k^p(p) of the Krawtchouk polynomials are 
independent of N. 

By a limiting argument from Theorem 16.31 we can obtain a closed formula for 
the connection coefficients in the case of a cyclic permutation. 
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Theorem 8.5. Let pi,... ,pd G (0,1) and \p\ = pi + ... + pa < For v, p G Nq 
and |u| = 1^1 = n, 

(8.13) 


where p = [pi,..., ’pd-i) with 
(8.14) p, = 


X . Vd-i)] p, n), 

PiPj+i 


(1 - Pi)(l + Pi — |Pj+il)(l + Pi - \pj\) 

Furthermore, we also have 

(8.15) 

{{pd,Pd-i ■ ■. ,/i2);p,n), 

where p = {pi,..., Pd-i) with 

/o ic-i n = _ PiPd+i-j(l ~ IpI) _ 

^ - ^ {i-\p^_^\){i-\p^^^_^\)ii + p,-\p\y 

Proof. Setting k = t{pi,..., pd, 1 — |p|) in equation (16.311 in Theorem 16.31 it is easy 
to check that 

l3o = tpi and A = A + Pi - \Pd+i-j\) + J = l,2, ...d. 

Taking the limit t —>■ oo in Rp(9; ft, n) and ignoring positive factors independent of 
V we obtain 

d-l 




i=i 


X 2 ^"! 


-pd+l-jAp'^ J+2|_|,y<i+l (1 - |Pd_j|)(l + Pl - |Pd+i_j|) 

,d-i+2| _ \ud-j\ 


Changing the product index from j to d — j, the above expression becomes 

UdJ* ■(1-I(.,+4)(1 + P1-|P,I)J' 

Applying now the Pfaff identity (see [TJ Theorem 2.2.5]) 

' a, b 


2 F 1 \ ;x 


a IP I ^ ^ 

= (1-X) .F,( ^ 


with a = -pj+i, b = |p^+^| - c = = -n+ + Wj-i\ and 

ignoring again positive factors independent of u we obtain 

d-l 

(8.17) (-l)l'^^l 


/U-j+l 


i=i 


X 27^1 


PiPj+i 


-Pj+i,-Vj . (1 ~ |Pjl)(l + Pi - |pj+il) 

^-n+ |p^'+2| + \vj_i[ 

The parameters pj defined in (18.141) satisfy 

1 - Ipj-iI ^ (1 - |Pjl)(i + pi - Ipj+iI) 

pj 


P 1 P 3+1 
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and comparing the last two formulas with (I8.3|) shows that, up to a positive fac¬ 
tor independent of v, equation ()8.17|1 coincides with the orthonormal polynomial 
• I ^d-i)]P, n). Furthermore, taking the limit in the weight func¬ 
tion wr{v,I 3, n) in (16.31) and using the formula (13.2L we obtain 

'^d^pT f}^ vj'- (i + pi - 

which can be simplified to 

{i-\p\Y- _^_ 

vY uY p'YY + pi- \Pj+i\Y^+''^+^' 

Up to a factor independent of v, this is precisely the weight function K^„(pi,..., Vd-i 
as can be directly verified using (18.141) and 


(i-pi)(i + pi - IpI)' 

The latter can be directly verified by writing pj in partial fraction. Positive multi¬ 
plicative factors independent of the variable v are not accounted for in the above 
computations, but they are unessential since an analogue of (I2.12|) for kj ^ completes 
the proof of (j8.13ll . Applying now Corollary 18.31 we obtain equation ()8.15p . □ 


9. Orthogonal polynomials on the unit ball and the unit sphere 

For the unit ball = {x : ||x|| < 1}, where ||a;|| denotes the Euclidean norm of 
X in the classical weight function is given by 

W^{x) = {l-\\xfY, P>-Y a.el'". 

A more general family of weight functions on the ball is defined by 

d 

One orthogonal basis for V^(1F®) is given by ■) ■ |q;| = n} with 


(9.1) P„"(VF®;x) := n(l - ||x,_if 


- + 5.«i + |) 


i=i 




where Xj = Xj{K,a) := -|- -I- d — j and denotes the generalized 

Gegenbauer polynomials that are orthogonal with respect to |a;p'^(l — on 

[—1,1]; see O p. 266]. It should be noted that the parameters in our weight function 
differ from those in by an additive 1/2. 

The orthogonal polynomials on the unit ball are known to be related to orthog¬ 
onal polynomials on the simplex. In fact, if {P/f : juj = n} is a basis of Vn{WY) 
on the simplex, it is known m Section 4.4] that {P^{xY ■ ■ ■, ■ \v\ = n} is a 

basis for the subspace of V)^(IF®) that consists of orthogonal polynomials that are 
invariant under sign changes. The proposition below completes the relation in the 
other direction. 
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Definition 9.1. Let : \v\ = n} he a basis o /on the simplex. Let e £ 
{0, l}'^ and e' = (e, 0). Assume {n — \e\)/2 G No- For u G Nq with \v\ = {n— \e\)/2, 
define 

Qt,eix) ■■= 

It is evident that is a polynomial of degree exactly n. To simplify the 
notation, we write below e instead of e', i.e. for e G {0,1}'^ when we need to work 
in we think of e as (e, 0). 

Proposition 9.2. If {P^f : \v\ = n} is an orthogonal basis afV^iW^), then the set 
Xn ■= {Que ■ W\ — S No,e G {0,1}'^} is an orthogonal basis for V^{W^). 

Proof. For n — |£| being an even integer and /3 G Ng with \(3\ < n, the integral 

Iil3)-=f Q'^^,ix)x^W^{x)dx= [ x^P;f~^%xl,...,xl)x^W^ix)dx 

Jm’‘ Jb‘‘ 

is equal to zero if /3 — £ contains an odd component, as can be seen by changing 
variable Xi 1—>■ —Xi. In the remaining case, we can write j3 = £ + 27 with jyj < , 

so that 

/(£ + 27 ) = f x^"<P!f+%xl, ..., x 1 )x'^‘^W^{x)dx. 

Since x^‘^W'^{x) = WK.+e{x \,..., x^^)xi ■ ■ ■ Xd, it follows from the identity 

[ f{xl,...,xl)dx= f f{xi,...,Xd)—=J^= 

that I{e + 27) = 0 by the orthogonality of since I7I < \v\. The mutual 

orthogonality of Q” ^ in follows from the parity of x^ . Finally, it is easy to check 
that the cardinality of is equal to dim V^, so that is a basis of V^(kF®). □ 


If P^{x) are the Jacobi polynomials on the simplex, defined in (12.71) of Section 
[21 the basis in the proposition above is exactly the basis (19.11) on the unit ball. 


Proposition 9.3. Let Pf; be the Jacobi polynomials on the simplex. Then, up to 
unessential constant factors, the set in Proposition \9.2\ coincides with the basis 
{P;^(IT»:-);H=n,aGNg}. 


Proof. The generalized Gegenbauer polynomials are related to the Jacobi 

polynomials by (Bp- 25-26]) 




(A + fj,)r 
id + ^)r 
(A -|- 
id + k) 


P^ =^2^2-1) 


m+1 


tP^ 


I). 


Using the above two identities, dSIl), dnn) and that Oj {k -|- £, u) = Xj {k, 2 v + e) for 
Qj dehned in (|2.6|) , one can verify that 


(9.2) 


"p;+^(a 


)^Pl^^{W^;x) 


where the symbol = means that the equation holds up to a multiple constant. □ 








CONNECTION COEFFICIENTS FOR ORTHOGONAL POLYNOMIALS 


27 


Let T be a permutation in Sd- Let tk = (t(ki, ..., Kd), Kd+i). For x € the 
weight function LF® satisfies Kd+i)('^ 2 /) = ^*( 2 /)- to see that the set 

{P”(fF®^; rx) : |q;| = n} is also an orthogonal basis of V^(fF®). Consequently, we 
can consider the connection coefficients 6 ^ ^(k) defined by 

|/3|=n 

From the above discussion, it is evident that bj, ^{k) can be derived from cj, ^(k). 

Theorem 9.4. Let e G {0,1}'^. Then, for v G Nq, the connection coefficients 
^ 2 u+e ~ ^ if P f £ for some fi gNq and 

(9-3) ^ 2 ’'!^+e. 2 ^ 1+6 (s^) = + e). 

Proof. That = 0 when f does not have the same parity as 2n + e follows 

directly from Proposition 19.31 using the connection coefficients of cf ^{k). The case 
(19.311 follows as a consequence of (19.211 since we can ignore the irrelevant constant 
factors when we consider normalized connection coefficients. □ 


In the above consideration we assume t G Sd instead of 5^+1 to exclude the 
permutation in 5’^+! acting on Kd+i, because such a permutation will have to 
act on Xd+i ■= y/l — ||a;|p as well in order to keep the weight function invariant. 
However, Xd+i is not a polynomial. To include permutations that act on Kd+i, we 
need to go back to the orthogonal polynomials on the simplex. Let us illustrate the 
situation in the particular case d = 2 for the classical weight function 

w«{x) = {i-xj-xir, /i>-i 

on the disk B^. There are two orthogonal bases that are well known O Sect. 2.3]. 
The first basis of V^(LF®) is (I9.1|l . in Cartesian coordinates, which we reindex by 


(9.4) P®Jxi,X2) = - xDic'; 


,/i+ 2 


X2 




0<j<n, 


whereas the second one is given, in polar coordinates {xi,X 2 ) = {r cos0 , r sin6) , by 


Qlfxi,X 2 ) = cos(n - 2 j) 0 , 0 < j < f, 

Q] 2 ixi,X 2 ) = sin(n - 2j)9, 0 < j < f. 

Recall the orthogonal polynomials P!f_j j, defined in (14.11) . with respect to the 
weight function Wk,{xi,X 2 ) = — xi — X 2 Y^ on the triangle T^. The first 

basis {P®„ : 0 < j < n}, as shown by Proposition 19.31 arises from Pfi-jj with 


Proposition 9.5. The basis (19.511 agrees, up to a multiple constant, with the basis 
in Provosition 1 .9. S\ that arises from : 0 < j < n} with k = (— i, — i, /i). 


Proof. For n = (—i, — i /r), we obtain by (14.41) that, for 0 < ^ < n, 


d-,’,7(4.1^) = (-i)"-pyjr'(2x; + 2x» - ixxf+x^)'p; 


\n—i 


2 \f ui-h-Y 


2x2 


X? + Xn 
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Since ^ (cos 0) is, up to a multiple constant, the Chebyshev polynomial 


Ti(cos0) = cos 10, and — 1 = 2sin^0 — 1 = — cos20, written in polar coordi¬ 

nate gives 

- I)r 2 ^cos 2 ^ 0 , 

where we use = to denote that the identity holds up to a constant, which is evidently 
Q^-i i{xi,X 2 )- Similarly, for k = (i, ^,/r), we obtain by (14.4|) that, for 0 < ^ < 
n — 1, 

xiX2pI^^I’!^-^^(^(xI,x\) = sin20P^^ = ’ = ^(cos20). 

Since P^'^'^\cqs20) = ^'^^^^ 2 ?^ ’ follows that the above polynomial is, up to a 
multiple constant, 2 (xi,X 2 )- This completes the proof for the basis of even 

degree. Using the identities 


, sin(2£-fl)0 (-ii), , cos(2£-fl)0 

PP' ='’(cos20) = — and P^ =’ = ^(cos20)= ^ ’ 


sin 0 

the proof for the basis of odd degree follows similarly. 


COS 0 


□ 


For the two bases in (19.41) and (19.51) , we can define the connection coefficients by 

n 

Ql-],i{Xl,X2) = '^bYmPZ,n{Xl,X2), 1 = 1,2, 

m —0 

where 0 < j,m < n and the hat symbol again denotes orthonormality. These 
coefficients can be deduced from the connection coefficients in Section |4l as shown 
by Theorem [01 For example, since P 2 m, 2 nixi,X 2 ) = Pn-m,mixi:xl) and Q'^^ i = 
Pn-e’e(xi,X 2 ), where k = (—^,—^,fi), it follows by the definition of cJ^(k), 


/52n _ 

Wn-e,l — X^n-i,e 


V 

/ . 4,m 

m—0 


{^)Pn-m,m ~ 


V 

m=0 


(k)P, 


2m,2n' 


Hence, 'bl'^ 2 m+i — ^ and 'bl'^ 2 m — particular, 'bl'^ 2 m ‘^an be written in 

terms of a Racah polynomial. The other cases can be worked out similarly. 

A similar result holds in higher dimension for the weight function W^. The 
definition of the corresponding basis in the polar coordinates will involve spherical 
harmonics dehned on the unit sphere of which are homogeneous polyno¬ 

mials orthogonal over with respect to the surface measure da. More generally, 
we can consider the weight function 

d+l 

Wk(x) = PJ Kj > -1, X & 

Evidently WK.(x)da becomes dcr if ki = ••• = Kd+i = — 5 . Let 'H‘^^(wk) be the 
space of homogeneous polynomials of degree n orthogonal with respect to Wk. on 
It is known that an orthogonal basis for can be derived from orthogonal 

polynomials on the unit ball. More precisely, for y G write y = r(x,Xd+i)', 

then a; € and (a;,a;d+i) G Let ip{x) = y'l — ||a;P. Define, for u G Nq, 
homogeneous polynomials 

yW( 2 /) |u| = n, 

■=r'^Xd+iP!;^~^((pW^;x), |u| = n - 1 . 


(9.6) 
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Then the set : \v\ = : \v\ = n—1}, restricted to is an orthogonal 

basis for (see, for example, O Section 4.2]). 

Comparing the above with Definition 19.11 we can add another parity at the last 
variable and establish the following result. 

Definition 9.6. Let {Pl^ : \v\ = n} be an orthogonal basis of V^(V14;) on the 
simplex. Let e £ {0,1}'^“'"^ and y = r{x,Xd+i). Assume (n — |e|)/2 £ No- For 
€ Nq with jj^l = (n — |e|)/2, define 

Ql.eiy) [xi,... ,Xd+i) e 

It is not difficult to see that Qf ^ is a homogeneous polynomial of degree n in 
y. Comparing with Definition 19.11 we see that the following proposition holds as a 
consequence of Proposition 19.21 and (19.61) . 


Proposition 9.7. The polynomial is a homogeneous polynomial in y and the 
collection of all such polynomials restricted to the sphere, {Qf,, '■ s £ { 0 ,1}'^^^, jz^j = 
(n — |e |)/2 £ No}, is an orthogonal basis of {w k). 


Note that, restricted to the sphere, Xd+i = \/l — ||a;|P is a polynomial of degree 
1. When Pff are the Jacobi polynomials on the simplex and (a;,a:d+i) is written in 
the standard spherical coordinates, the basis in the above proposition is precisely 
the one given in [5l Theorem 7.5.2]. In particular, this gives the standard basis for 
the spherical harmonics when ki = ... = Kd+i = — 5 , see, O Theorem 4.1.4]. 

As an example, let us state explicitly this basis for Ti^ of spherical harmonics 
of three variables in terms of the orthogonal polynomials on the simplex. If Y 
is a homogeneous polynomial, then Y{x) = r^Y(x') for x = rx' with r > 0 and 
j|x'j| = 1. Hence, we only need to give the basis in terms of x'. 


Example 9.8. Orthogonal basis for TL^ of the spherical harmonics on S^. Let 
: 0 < fc < nj be the Jacobi polynomials defined in (14.11) . Then, 
for {xi,X 2 ,X 3 ) £ S^, 

j <n}U {xiX 2 Pit’fl~]\xl,xl ),0 < j < n - 1} 

U {xiX3Plt\Jf]\xl, xl), 0 < } < n - 1} U {x2X3P^Zliyjf {xl, xl), 0<j<n-l} 
is an orthogonal basis for TL^ni 

{xiPj,l'~/'~^\xl,xl),0 < j <n}U {x 2 PlZ]f~^\xl,xl), 0 < j < n} 

U {x3P^Z]^'~^'^^ (xi, xl), 0 < J < nj U {xiX2X3P^Z'fZ]jixl, xl), 0 < } < n - 1} 
is an orthogonal basis for the space TLln+i ■ 


To see directly that this basis is indeed the classical basis of spherical harmonics 
on S^, we need to use the standard spherical coordinates and carry out computations 
as in the proof of Proposition 19.51 We leave it to the interested readers. 

The restriction of {QlZ,{T{x,Xd+i)) : e £ {0, l}'^'*'^, (n— l£l )/2 G No} is evidently 
another orthogonal basis of TLfiWfZ)- We can define the connection coefficients 
between this basis and the basis {Q^g}- There is no more mystery or work needed 
to be done in this case. These connection coefficients can directly be computed 
from those for the Jacobi polynomials on the simplex. 
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